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Abstract 

We calculate the scattering intensity of dilute and semi-dilute solutions of 
star polymers. The star conformation is described by a model introduced 
by Daoud and Cotton. In this model, a single star is regarded as a spheri- 
cal region of a semi-dilute polymer solution with a local, position dependent 
screening length. For high enough concentrations, the outer sections of the 
arms overlap and build a semi-dilute solution (a sea of blobs) where the inner 
parts of the actual stars are embedded. The scattering function is evalu- 
ated following a method introduced by Auvray and de Gennes. In the dilute 
regime there are three regions in the scattering function: the Guinier region 



(low wave vectors, qR <^ 1) from where the radius of the star can be ex- 
tracted; the intermediate region (K gii: < that carries the signature 
of the form factor of a star with / arms: I{q) ~ q~^^/^; and a high wavevector 
zone {qR » /^''^) where the local swollen structure of the polymers gives rise 
to the usual decay. In the semi-dilute regime the different stars interact 
strongly, and the scattered intensity acquires two new features: a liquid peak 
that develops at a reciprocal position corresponding to the star-star distances; 
and a new large wavevector contribution of the form q~^/^ originating from 

the sea of blobs. 
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A star polymer is a branched macromolecule with a small central core from where emerge 
several end-attached linear chains |jl|] . Not only can star polymers be regarded as prototypes 
of branched systems, but also can solutions of star polymers be considered as model systems 
to study polymer colloidal stabilization or the effects of internal segment density on the 
behaviour of viscoelastic solutions. Experimental work in dilute solutions includes studies 
of viscoelastic behaviour and chain conformation by rheology , light scattering and 
small angle neutron scattering (SANS) 0. The chain dimensions and structure factors in 
semidilute solutions or melts has also been studied by SANS • Theoretical work on the 
conformation of star polymers in good solvents was pioneered by Daoud and Cotton and 
Birshtein and Zhulina Renormalization group calculations |12|-[T^, Monte Carlo |1T5|JT6 



and molecular dynamics simulations |lT7| , p!8| corroborate the Daoud-Cotton (D-C) picture 
at least in the dilute regime where the simulations have been performed. However, only a 
limited amount of information is available to describe the scattering intensity /(q) of star 
polymers. Analysis of the scattering data in the litterature is performed with postulated 
expressions or limited to the Guinier regime from where the star radius of gyration can be 



extracted W^ . In this paper we study the different factors which contribute to the scattering 



intensity of a star solution at arbitrary concentrations. We follow a method introduced by 
Auvray and de Gennes pO| , pl| to evaluate the intensity scattered by an adsorbed polymer 



layer. This method allows for an explicit computation of the scattered intensity of a star 
solution in the full wavevector range < g < {b{b^ / v))^^ , where b is the Kuhn length and v 
the excluded volume parameter. 

The structure of a star polymer with / arms of polymerization index can be described 
by the Daoud Cotton |jlO[ model. Attachment of the chains to a central core effectively forces 



the local polymer density to be everywhere inside the star above overlapping concentration. 
The star can therefore be described as a semi-dilute solution p2| , p3| , with a local, position 
dependent screening length ^(r), where r is the distance from the center. Pictorially, we 
represent this by associating with each arm a string of blobs of increasing size ^(r). The 
radial dependence of the blob size ^(r) can be obtained by noticing that at a distance r 



from the center there are / blobs of cross section ^(r)^ occupying a total area of 47rr^. 
The blob size thus varies as ^(r) ~ rf-^/^ and the corresponding local polymer volume 
fraction as (psi^) ~ p^^ih/rY^^ . Note that there is a crowded region of size Rc — bf ^^'^ 
in the middle of the star where the concentration reaches one. The size of the star can 
be obtained from monomer conservation Nfb^ — An J^r'^ dr4>s{r). Neglecting the small 
core region one gets R ~ aN^^^p^^. Figure 1 sketches the structure of the star solution 
above the overlapping concentration 0* = fNb^/R^. The inner part of each individual 
star is still described by a position dependent correlation length ^(r) and the corresponding 
concentration profile (psi^) — However, the profile doesn't vanish at a distance 

R but levels off to a constant value 0p at a distance Rgd, determined from the condition 
(f)s{Rsd) = <Pp- The concentration (f)p is the polymer concentration in the region known as 
the sea of blobs, where the outer sections of the arms overlap. The stars embedded in the 
sea of blobs occupy a volume fraction $s and the sea of blobs the remaining fraction 1 — $s. 
Above the overlapping concentration 0*, the values of the polymer concentration in the sea 
of blobs 4>p, of the fraction occupied by the stars $s and of the radius of the embedded stars 
Rsd can be calculated from mass conservation: 



-1 5/4 



0p ~ and Rsd — R 



n3/4 



(1) 



Notice that an increase of the polymer concentration leads to an increase of the fraction 
of space occupied by the sea of blobs and a shrinkage of the inner star dimensions. 
The scattering intensity due to the chains in solution is given by 

/(q) =< a*(q)a(q) >, (2) 

where a(q) is the amplitude scattered by the J\f monomers in the system and <> denotes 
an average taken over all possible configurations. The scattered amplitude is written as 

= Yli'^i ~ expO'q.rj} (3) 

i 

with rii — ns the difference between the scattering lengths of the i*'*-monomer and the solvent. 



We assume hereafter that all monomers have equal scattering lengths, hq. Defining 0(r) as 
the monomer volume fraction at position r, expression can be written as 

-^(q) ^ Jf J d^rexp{jq.r} J d^r' < p(r')p(r + r') >, (4) 

with p(r) = {riQ — Us) (f){r) . The scattering amplitude p(r) of the semidilute solution of 
polymer stars described above can be represented as 



p(r) = E '^(r - i-i) ® [Ps(r)n(r)] + pp(r) 



1=1 



1 - ^(5(r - Ti) ® n(r) 



i=l 



(5) 



where Ng = M / {N f) is the number of micelles in the volume V, the symbol ® denotes the 
convolution integral /(r) ® gij) = J d^r'f{r — r')g{r') and n(r) the step function n(r) = 
1 for r < i? ; n(r) = for r > i?. The first term on the r.h.s. of equation ^ effectively 
locates Ns spherical inner parts of stars at positions rj. The second term accounts for the 
monomers in the sea of blobs. The scattering amplitudes of the inner parts of the star Ps(r) 
and of the sea of blobs Pp(r) can be expressed as the sum of an average contribution and of 
a fluctuating part 

Ps{p){r) = {no - ns)(l)s{p){r) + {uq - n^)5</)s(p)(r) (6) 

By inserting expressions (^) and @ into equation (^) we evaluate the total scattering in- 
tensity under the following approximations: i) density fluctuations in different stars are 
uncorrelated. ii) density fluctuations in the sea of blobs are uncorrelated to density fluctua- 
tions in the inner part of the stars. This is a reasonable assumption except perhaps for the 
correlations between the outer blob layer of the stars and the neighbouring blob layer in the 
sea of blobs. The total scattering intensity can be written as a sum of three terms 

^(q) = {no - n,yV [/pp(q) + /s.(q) + /sp(q)] (7) 
The flrst term Ipp accounts for density fluctuations in the sea of blobs 

/pp(q) = (1 - < <^0j(q) >, (8) 
the second term Igs accounts for density fluctuations in the inner part of the stars 



/,,(q) = < 5(j)'M) >, (9) 

and the third term Isp expresses correlations between the different stars: 

/,,(q) = <l>,5(q)A0(q)2. (10) 

This last expression contains the stars structure 

factor S'(q) = l/NsY.m.n=i^w{j^i''^m — rn)} that expresses correlations between the po- 
sitions of the stars in the solution. It is weighted by the form factor of the inner parts of 
the stars A0(q)^ = (0p(q) — 0s(q))^- We do not attempt here a calculation of the structure 
factor S{q), our results provide instead a systematic approach for extracting ^(q) from the 
experimental data. 

The results presented above can also be used for dilute solutions of polymer stars. For- 
mally, one sets the sea of blobs contribution equal to zero and takes the limit S'(q) 1. 
The volume fraction occupied by the stars is in this limit $s = 0/0*. 

We now calculate explicitly the polymer contributions to the three different terms @- 
(p!0|). The contribution from polymer-polymer correlations in the sea of blobs can be written 
as the fourier transform of g{r), the pair correlation function for density fluctuations. This 
is a function which decays algebraically as r~^^^ and vanishes rapidly for distances larger 



than the correlation length = ^{(pp): g{r) ~ exp{— r/(-i/27/20^p)}/r'^/^. The numerical 



value y 20/27 is chosen such that the correlation length can be directly measured from the 
lorentzian decay of the structure facture at low wave vector < 50p(q) >~ (1 + (o'^p)^)^^- 
Performing the fourier transform of the correlation function leads to 



< (50j(X) >= Cte 



b 



nl/3 



(5 sin {I arctan y^27/20X} 



(11) 



V 3 (X3 + 27/20X5)1/3 
with the dimesionless wavevector X = q^p. For small wavevectors {X ^ 1) one has by 
construction < 50p(X) >~ Cte(^p/6)^/^(l + X'^)~^ while for large wavevectors (X ^ 1) one 



gets < 502(X) >~ Cte^5/4 {2<d/2lfl\ip/b)-^l\'''l\ 

The contribution from the density fluctuations inside the stars can be regarded as the 
scattering from a semidilute solution with a position dependent correlation length: 



< 64>1{X) >= Cte 



b 



1/3 



/5 ri sm{^arcta^^^^_ 
3 Ja^ ^ {{Xy)3 + 27/20(X|/)5)i/3 



27/20X?/} 



(12) 



with a = Rc/Rsd, the ratio between the compact core radius Rc and the inner star radius 
Rsd- The form of this contribution is similar to the contribution from the sea of blobs. For 
instance, in the large wavevector limit {X ^ 1) one has < 50^(X) >= 9/5 < 6(f)p{X) >. 
The form factor A0(q)^ in equation ([T0|) is calculated from 

1 2 



A0(q) 
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2 , smgr 
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(0(r) - (f)p 



(13) 



with (/)(r) = 1 for r < Rc ; (l){r) = /^/^(fe/r)"^/^ for Rc < r < Rsd and (pp a constant over all 
the integration range. Performing the integral leads to 



'h' 


8/3 - 
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a^/^,(aZ)-,(Z) + ^£, 



dy 



smy 



(14) 



where we defined a dimensionless wavevector Z = qRsd- The asymptotic behaviour of the 
form factor for large wavevectors {Z ^ 1) is A0(Z)^ ~ R^J^^C,p^^^q^^^^^ ■ 

Comparison of the three factors described above show that the scattered intensity is 
dominated at high q vectors by the q~^/^ behaviour that is due, in dilute star solutions, to 
concentration fluctuations inside the star and, in semidilute star solutions, to concentration 
fluctuations in both the inner part of the stars and in the sea of blobs. At low wave vectors 
we predict for dilute solutions a Guinier regime followed by a q^^^/^ slope in a log-log 
representation, due to the form factor of the stars. For semidilute solutions, correlations 
between different stars are important, and the form factor multiplies the structure factor 
which measures those correlations. In practice this must lead to the appearence of one or 
more correlation pics in the Guinier zone or above. The crossover between the q~^/^ and the 
g-io/3 j-egimes is expected to occur at wave vectors qRsd ~ f^^^ in semidilute solutions and 
at qR ~ /^/^ in dilute solutions. 

The general form of the scattering intensity of a polymer star solution, as outlined above, 
is in good agreement with results from both experiments and simulations. In the following 
paragraph we quantitatively confront our predictions to experimental data from reference [§. 



The stars considered in reference are synthesized by the so-called core-first method. In 
this method, a core of divynyilbenzene (DVB) is synthesized and from the surface of the 
core, anionic polymerization of polystyrene takes place. This method has the advantage of 
providing high functionality objects, but with the inconvenient of a high polydispersity of 
the functionality (polydispersity of core masses). Also, since the linear chains are grown 
from, sometimes, massive DVB nodules, steric repulsion between monomers in the early 
stage of the polymerization are probably responsible for branches which are longer than the 
value targeted during synthesys. In the present case, the branch mass was targeted to be 
M„ ~ 2 X 10^. 

We plot in figure 2 a typical scattering curve from a semidilute star solution. The inter- 
mediate and large q ranges clearly exhibit slopes —10/3 and —5/3 as predicted. In order to 
extract the structure factor of the stars solution we first subtract the large wavevector con- 
tributions. The unknown constant of equations (^TJ) and (|1^) can be obtained from data on 
semidilute solutions of linear polymers. Indeed, the scattering at small angles of linear poly- 
mers in a semi-dilute regime is given by an Orstein-Zernicke equation containing the same 
multiplicative constante of equation (^) and (|T2|) . In our case we extracted the numerical 
value of the constant from data on polystyrene in deuterated toluene We fixed the ratio 
a at the value 0.05. The curves obtained after subtraction are presented in figure 3. For the 
more diluted sample, = 0.055, the slope of data in a log-log representation approaches the 
theoretical slope of —10/3. Deviations from this slope for the more concentrated samples 
migth have several origins. Second or higher order pics in the structure factor would fall in 
this range, therefore introducing oscillant deviations from a pure linear slope in the log-log 
plot. Also, the quality of the stars may be here of importance, particularly the size, shape 
and functionality of the star cores. A few large cores with a small number of attached arms 
would also contribute significantly to the scattered intensity in the intermediate range. 

As a final step we divide data from figure 3 by the calculated form factor of equation 
(|I^). The structure factors obtained by this method are plotted in figure 4 for five different 
polymer concentrations. The functions obtained display the usual features of liquid structure 
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factors, in particular the position of the peak increases with concentration as (j)^/'^ [|]. It is 
also important to stress that the procedure of subtraction is extremely sensitive to both the 
chemical quality of the stars and the correspondent quality of the scattered data. We expect 
that structure factors of better quality can be extracted from data of a carefully defined 
chemical nature. Such structure factors could then be confronted to existing theories of 
liquid structuration. 

In this paper we have calculated the scattering intensity of a star-polymer solution, both 
in the dilute and semidilute regimes. The structure of the solution and the conformation of 
the stars has been described by the Daoud-Cotton [O] model, and the intensity calculated 



following a method introduced by Auvray and de Gennes |]20| , |21| . Our results show that 
the scattering intensity can be described as the sum of three terms. The first term is the 
product of the structure factor of the solution by the form factor of the star. We did not 
evaluate the structure factor for a star solution, but our results provide a procedure that 
allows its extraction from experimental data. The form factor is the fourier transform of 
the average concentration profile of the star. This term gives rise to a caracteristic q^^^/^ 
wavevector dependence in the intermediate q range (1 ^ qRsd ^ P^^)- The second and 
third terms are independent of the structure factor and measure concentration fiuctuations 
in the inner parts of the stars and in the sea of blobs. They have both a q~^f'^ dependence, 
and dominate the scattering spectra at high wavevectors qRsd ^ ■ Our data is consistent 
with available experimental data on star solutions. 
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FIGURES 

FIG. 1. Schematic representation of a semi-dilute solution of star-shaped polymers. The aver- 
age concentration profiles between two neighbouring stars is also represented. 

FIG. 2. Scattering intensity as a function of the scattering vectror in a log-log scale for a solution 
of star polymers. The concentration is 50mg/cm^ and the average functionality (overestimated) is 
130. The star branch mass targeted in synthesis was M„ = 2 x 10^. Straight lines of slopes —5/3 
and —10/3 are also represented in the figure. Data from ref [8]. 

FIG. 3. Scattering intensity as a function of the scattering vector in a log-log scale for solu- 
tions of very polidisperse star polymers at different concentrations when q~^l^ terms have been 
subtracted. The slope —10/3 is also shown for every data set. 

FIG. 4. Structure factors of solutions of polydisperse stars at different concentrations obtained 
from the procedure described in the text. 
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